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SEASONAL CHANGES IN TI1E COMMUNITY STRUCTURE OF THE 
MACROBENTHOS INHABITING THE INTERTIDAL SAND 
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1 he University of Chicago, Committee on Evolutionary Biology, Chicago, Illinois; 
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Since the pioneering works of Peterson (1918) and Thorson (1957), most 
of the attention of workers in the field of marine benthic ecology has been di¬ 
rected to short-term sampling surveys of particular environments in an attempt 
to delimit given macrobenthic species assemblages or communities. These investi¬ 
gations have generally concluded that given species assemblages are greatly in¬ 
fluenced by sedimentary parameters ( e.g ., grain size), and that the communities 
are not tightly functioning units but rather are composed of individualistically dis¬ 
tributed species which respond to a complex set of environmental and biological 
parameters. The result of such findings is that given particular sedimentary 
parameters, a posteriori predictions of what faunal assemblages will occur at a 
particular site are possible. However, over longer periods of time, the predictive 
value of determining the distribution of given species groupings is greatly reduced. 
Any naturally occurring community of organisms is not only distributed in space 
but also in time. To better understand patterns of community organization, there¬ 
fore, more long-term studies are needed to provide basic information regarding 
temporal community structure and to elucidate those factors affecting species 
compositional change. 

Studies on temporal changes in community structure and species composition 
are not common. Of these studies (e.g., Sanders, 1960; Pearson, 1971, Tenore, 
1972; Lie and Evans, 1973; Watling, 1975; Levings, 1975), most are concerned 
with subtidal environments. Surprisingly, there are fewer studies dealing with 
time-related community changes in intertidal environments ( e.g Bloom, Simon 
and Hunter, 1972; Johnson, 1970; Holland and Polgar, 1976). Intertidal environ¬ 
ments provide an excellent opportunity to study the effects of seasonal change, 
since physical variables fluctuate with greater amplitude than in bordering subtidal 
habitats, thus subjecting organisms inhabiting them to greater amounts of stress. 

The purpose of this study is to describe the general macrobenthic community 
structure of a portion of the intertidal sand and mud fiats of Barnstable Harbor, 
Massachusetts, and also to examine in greater detail seasonal changes in patterns 
of community structure and species composition at one of the sedimentary environ¬ 
ments. This study is presented to provide a context for a more detailed examina¬ 
tion of methods of resource allocation in deposit-feeding faunal assemblages and 
micro-dispersion patterns of surface-feeding polychaete annelids. 

1 Present address: Department of Biology, Woods Hole Oceanographic Institution, 
Woods Hole, Massachusetts 02543. 
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Materials and Methods 

Barnstable Harbor, located on the northern shore of Cape Cod, is a shallow 
tidal enbayment about 5.5 km long and 1.2 km wide. Over the past several 
decades a relatively large amount of research has been directed towards various 
aspects of the marsh ecology (see Redfield, 1972, for a partial list of references 
and an excellent description of the general development of the marsh system). 
The mean tidal range is about 2.9 m and approximately 60% of the tidal flats, 
where this study took place, are exposed at low tide (Ayers, 1959). 

The physical environment of the region is typical of intertidal regions along the 
southern New England coast of North America. Water temperatures range 
from about 20° C during July to September to about 0-1° C during December 
to February. Sediment temperatures range from about 30° C in August to 1-2° C 
in December to February. The winter of 1974-1975 was unusually mild, and 
accumulation of ice on the flats was never observed. Ice buildup on the flats 
normally occurs from January to February (Green and Hobson, 1970; personal 
observations in the winter of 1975-1976). Ayers (1959) reports that salinity near 
the present study area varies between 29.6 to 31.7%, with very little freshwater 
drainage in the harbor. 



Figure 1. Map of western Barnstable Harbor, Massachusetts, indicating sampling localities. 
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The currents of Barnstable Harbor are controlled almost exclusively by the 
tides. Due to the large tidal range, strong currents occur in the harbor. Both 
Ayers (1959) and Sanders, Goudsmit, Mills, and Hampson (1962) have 
measured the currents at different regions within the harbor and found maximums 
from 30—80 em/sec occurring during flood tides. The strongest currents occur 
just before the time of maximum ebb tide. The flats within the harbor, however, 
are not greatly affected by the strong currents, since they are always exposed 
before the time of maximum ebb current. 

Changes in sediment accretion and erosion were monitored over a yearly period 
at stations 1 and 2 (Fig. 1). Results indicate minimal changes in the level of 
'substratum at these localities during a 15 month period (on the order of 1-3 cm). 
Maximum sediment fluctuation occurred in the spring. Larger amounts of sedi¬ 
ment movement were common, however, in the tidal creeks and areas closer to 
the mouth of the harbor. 

Sixteen preliminary stations were sampled in the western portion of the 
intertidal sand and mud flats in June-July 1974. Faunal data collected from these 
stations were used to help in delimiting the general characteristics of the various 
macrofaunal assemblages of this region. One station was chosen from the initial 
sampling program for more detailed long-term sampling. 

Samples were collected at low tide with a small hand held core 12.5 cm in 
diameter. The core was gently pushed into the substrate to a depth of 17-20 cm 
and removed with a shovel. The location of sampling sites at each station was 
determined by random sampling coordinates. Samples were brought to the labora¬ 
tory and gently washed through a 250 /mi mesh screen. The residues were pre¬ 
served in 70% alcohol, sorted under a dissecting microscope and the organisms 
encountered were identified to species. 

A small core (3.5 cm diameter) was used to sample surface sediment for sub¬ 
strate analysis. Sediment size-distributions were determined for each sample 
using standard wet-sieving analysis and mean particle size, and sorting coefficients 
were graphically determined from formulae of Folk (1966). The percentage of 
oxidizable organic carbon was estimated by the Walkley-Black wet oxidation method 
(Morgans, 1956). 

Faunal data were subjected to classification and cluster analysis in an attempt 
to delimit community structure, species composition and seasonal community 
dynamics of the benthic macrofauna. This type of analysis has received exten¬ 
sive attention in marine community studies ( e.g ., Lie, and Kelley, 1970; Stephenson 
and Lance, 1970; Hughes and Thomas, 1971 ; Stephenson, Williams and Cook, 1972 ; 
Boesch, 1973), and these authors have outlined the various methods. 

In this study, the per cent similarity coefficient (a derivative of the Czekanowski 
similarity coefficient) was used. Stephenson ct al. (1972) have briefly reviewed 
the characteristics of the different classification measures. The per cent similarity 
measure tends to be biased toward the more dominant and ubiquitous species, so 
that rarer species contribute little to the index. Preliminary data analysis com¬ 
paring the fanual data with and without dominant forms, however, revealed that 
the clusters were not appreciably altered. 

Two clustering strategies were initially employed: first, the more common 
group-average or unweighted pair-group method; and secondly, the newer and 
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more popular flexible sorting strategy. Lance and Y\ illiams (1967) and Williams 
(1971) have reviewed the basic properties of the two strategies. After initial 
data processing, it was found that the group-average clustering technique yielded 
the more instructive classification. The “flexible” method is a strong sorting 
strategy (William, 1971) which increases the distance between groups as the num¬ 
ber of elements which they contain increases. While this heightens differences 
between groups, it tends to divide these elements into groups that may obscure the 
important characteristics of individual systems. (Comparisons of the two clustering 
strategies are available upon request.) 

Both normal and inverse classification analyses were used in the study. Normal 
analysis (sometimes called Q-analysis) treats samples as individuals, each being 
composed of a number of attributes—the various species from a given sample. 
Normal analysis is instructive in helping to ascertain community structure and to 
infer specific ecologic conditions between sampling stations from the relative dis¬ 
tributions of the species. Inverse classification (termed R-analysis) is based on 
species as individuals, each of which is characterized by its relative abundance in 
the various samples. This type of analysis is commonly used to reveal species 
groups within particular habitats. 

Shannon’s formula (Shannon and Weaver, 1963) expressed as: IT = —2pi 
logopi, where pi is the proportion of individuals belonging to the i th species, was 
used to estimate species diversity. Evenness, the distribution of the individuals 
among the various species, was computed as: J' = H'/logo s, where s is the number 
of species in the sample (Pielou, 1966). 

Results 


Sediments 

A summary of the sedimentary features of the sampling localities is given 
in Table I. Four distinctive substrate types are evident: sandy sediment (stations 
10-16), muddy-sand sediment (stations 1-2), muddy sediment (stations 3 and 
5-9), and gravel-sand sediment (station 4). 

The muddy stations are all located in slightly higher and more protected areas 
of the marsh and are relatively free from tidal motion. These sediments are char¬ 
acterized by both high silt-clay fractions and organic carbon content. Dead and 
decaying fragments of Spcirtina alterniflora are commonly found throughout these 
sediments. 

The sandy stations can be further subdivided into two categories: fine (sta¬ 
tions 12-14) and coarse (stations 10, 11, 15, 16) sand. Coarser sand is generally 
found closer to the tidal channels. The sands in these areas are moved freely 
by tidal currents, and the surface is commonly marked by sand ripples indicating 
sediment instability. The finer sand stations are further from the channels 
where current velocities are lower. These sediments are generally well-sorted 
and contain low percentages of organic carbon. 

The very large sand flats north of New and Sand Islands (Fig. 1) are com¬ 
posed of well-sorted muddy-sands. These areas are more protected from current 
flow across the flats and any surficial ripple marks are normally the result of 
wave action. These sandy flats have slightly higher amounts of organic carbon 
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Table I 


Sedimentary and trophic group parameters from the sampling stations at Barnstable Harbor (see Figure 
1 for station locations). 


Station 

number 

Per cent 
silt-clay 

Mean phi 

Phi of 
sorting 

Per cent 
carbon 

Deposit-feeders* 

Suspension-feeders* 

Number 

species 

Number/ 
m 2 ( X 103) 

Number 

species 

Number/ 
m 2 (X 103) 

1 

11.9 

2.48 

0.18 

0.48 

19 

45.7 

2 

1.30 

2 

20.9 

2.47 

0.30 

0.78 

12 

39.0 

0 

— 

3 

91.0 

4.09 

0.40 

1.52 

3 

52.2 

0 

— 

4 

1.3 

0.53 

0.83 

0.37 

12 

6.3 

1 

0.08 

5 

48.6 

3.10 

0.33 

0.57 

15 

33.3 

0 

— 

6 

34.2 

3.00 

0.33 

0.66 

9 

18.7 

0 

— 

7 

88.7 

4.22 

0.45 

1.49 

6 

35.8 

0 

— 

8 

87.6 

4.10 

0.38 

1.54 

11 

21.1 

0 

—. 

9 

80.4 

3.81 

0.37 

1.35 

11 

23.5 

0 

— 

10 

4.1 

1.90 

0.35 

0.10 

8 

2.4 

3 

2.00 

11 

2.5 

1.78 

0.40 

0.07 

4 

0.9 

1 

0.20 

12 

11.4 

2.15 

0.48 

0.32 

17 

30.3 

2 

104.90 

13 

7.9 

2.13 

0.45 

0.33 

9 

11.9 

3 

44.80 

14 

7.5 

2.05 

0.43 

0.31 

9 

4.5 

2 

192.90 

15 

2.1 

1.47 

0.38 

0.22 

7 

6.5 

2 

28.20 

16 

2.6 

1.60 

0.32 

0.09 

13 

8.2 

1 

0.08 


* Feeding types determined from examination of gut contents, Sanders et al. (1962) and Barnes 
(1964). 

and silt-clay fractions when compared to the sand areas. Monthly sediment samples 
at station 1 indicated no significant temporal change in mean grain size or sorting 
coefficient according to analysis of variance (P > 0.05). 

The one station located on the south shore of the harbor consists of gravel, 
sand, and clay, all thought to he derived from erosion of the upland (Redfield, 
1972). This sediment is very poorly sorted, with low organic carbon content. 
The presence of patches of scattered stones and gravel characterizes this very 
heterogeneous substrate. 


General community description and species composition 

Table II lists the macrofaunal species occurring in 64 core samples obtained 
from the preliminary sampling of 16 stations in Barnstable Harbor. Of the 
47 species found in the survey, the majority were polvchaetes (29 species), 
crustaceans (9 species) and molluscs (6 species). Individual replicate samples 
taken at each locality were not lumped together for use in the classification analysis 
so as to obtain an idea of within-habitat sampling variability. The results of 
the inverse classification analysis used to delimit species groups, however, are based 
only on those species which occurred more than five times in the survey. Species 
eliminated from the analysis include ten which occurred once, five which occurred 
twice, and only one which occurred three times. 

Normal analysis. Figure 2 indicates that the samples are clustered into five 
major groups. Group A is composed of stations 12-15 containing clean, well- 
sorted sands with a low silt-clay fraction (2-5%). Group B, composed of samples 
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'Fable II 

Macrobeuthic species sampled at Barnstable Harbor , June—July t 1974 . 






Station 






Identifica- 

Species 1 









tion number 

1 

2 

.5 : 

4 

5 

6 

7 

jl| 

9 

10 

11 

12 

13 

14 

1 

15 

16 


Polvchaeta 

t 

strcblospio benedicti 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 


X 

X 

X 


X 

2 

Eleone heteropoda 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

3 

Lumbrinereis tenuis 

X 

X 


x 

X 


X 

X 

X 



X 





4 

Capitella capitata 

X 


X 

X 1 


X 


X 

X 

X 


X 


X 


X 

5 

Glyeera dibranchiata 









X 








0 

Polvdora ligni 



X 

X 

X 












7 

Phyllodoce mucosa 

X 

X 







X 



X 





8 

Ileteromastus fill form is 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 


X 


X 

<) 

Scolecole pules viridis 



X 


X 












10 

Scoloplos acutus 

X 




X 






X 

X 

X 

X 

X 

X 

11 

Spio setosa 

X 

X 


x 1 

X 

X 







X 




12 

Scoloplos robust us 






X 

X 










1.5 

Clymenella tor quota 





X ' 



X 









14 

Xereis nrenaceodonta 

X 

X 



x 



X 

X 



X 




X 

15 

Pectinaria gouldti 

X 

X 


X 

X 

X 











16 

Scolelepis squamata 

X 

X 










X 




X 

17 

Tharyx sp. 

X 

X 


x 








X 




X 

IX 

Xcphtys picta 

X 
















10 

Syllides verrilli 

X 

X 



X 

X 

X 

X 

X 

X 


X 

X 

X 


X 

20 

Spiophones bombyx 

X 

X 








X 

X 

X 


X 


X 

21 

Prionospio heterobranchia 

X 



X 





X 



X 





22 

Paraonis Jul gens 
















X 

2.5 

A r ic id ea je fl reys i i 

X 
















24 

Eleone longa 













X 


X 


25 

Xereis wens 




X 





X 








26 

Xepltlys caeca 










X 






X 

27 

A mphitrite ornata 




X 












28 

Pygospio elegans 









X 



X 

X 

X 


29 

Chaelozone sp. 














X 


Mollusca 
















.50 

Hydrobia sp. 

X 









X 

X 

X 

X 


.51 

V as sari us obsolctus 

X 

X 

X 



X 




X 





.52 

Gemma gemma 

X 

1 






X 


X 

X 

X 

X 


.5.5 

Mya arenaria 

X 

X 






X 

X 


X 

X 


X 

.54 

Relasa canaliculata 

X 

x 














.55 

Polynices duplicatus 

X 















C rustacea 

















.56 

Diastylis polita 















X 


.57 

Am pelt sea abdita 













X 




.58 

Oxyurostylis smithi 

X 
















.59 

Edothea triloba 

X 




X 



X 









40 

Crangon septemspionosus 








X 

X 








41 

Coroph i ti m i ns idiosu m 












X 





42 

Gant mar us sp. 










X 


X 





4.5 

Unicola irrorala 


X 















44 

Ampithoe rubricata 





X 

X 


X 









Hemichordata 

















45 

Saccoglossus kowalewskii 
















X 

Oligochaeta 


















46 

Phallodrilus monos per ma¬ 


















th ecus 












X 

X 


X 


47 

Monophylephorus irroratus 




.X 














collected from stations 10, 11, and 16, is also a clean sand environment. These sedi¬ 
ments, however, are coarser-grained with a smaller silt-clay fraction (2-3%). 
Group C contains the muddy sediment localities (stations 3 and 5-9) of Jules 
Island and Huckins Island. The major sedimentary features of these stations 
include high quantities of silt-clay (10-70%), small median grain size and high 
organic carbon content. The muddy-sand stations (1 and 2) fall into Group 
D. While sation 1 is more characteristic of a sandy sediment type, the larger 
silt-clay fraction makes it more similar to station 2 in sedimentary composition. 
Finally, Group E contains the very poorly sorted gravely-mud station. 

Inverse analysis . The 32 species found at the various localities in Barnstable 
Harbor were divided into 10 groups (Fig. 3) by the inverse classification analysis. 
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Groups 1, 3, 5, and 7 are comprised of species which characteristically occur in clean 
sand areas. While normal classification analysis divided the sandy stations into 
two separate groups, fine and coarse sand, species groups generated by the inverse 
analysis indicate that degree of within-habitat species fidelity is low. About 
39% of the species in these groups are found in both coarse and fine sands—the 
most striking visible difference being the presence of very large numbers of the 
small venerid bivalve Gemma gemma in the fine sands. 

Group 4 is composed of the most ubiquitous and numerically most important 
species of Barnstable Harbor. Streblospio benedicti, Eteone heteropoda, and 
Hetcromastus filiformis occur at each sediment type sampled. At sandy-mud and 
mud stations, S. benedicti and E. heteropoda were always ranked as the first 
and second most dominant species, respectively. 

Groups 8 and 10 consist of species normally associated with sandy-mud and 
mud stations. Although normal classification analysis separated these two sedi¬ 
ment types into discrete groups, the inverse analysis showed that between 
habitat species co-occurrence is high. Only two species, Scoleocolepidcs viridis 
and Scoloplos robust ns, in these two groups were restricted to just one sediment 
type. 

Group 9 consists of two species which are largely restricted to station 4, the 
most hetrogeneous substrate type. While the large nereid polychaete, Nereis 
virens, was also found at station 9, its greatest abundances occurred at station 4. 
The oligochaete, Monophylephorns irroratus, was never found to occur at any 
of the other 15 stations sampled. 

Groups 2 and 6 contain just one species each. Ampithoe rubricata, a small 
amphipod, and the epifaunal tectibranch, Retusa canalicidata, were found in the 
muddy-sand stations 1 and 2. Both species were relatively rare in abundance. 

In summary, inverse classification analysis produced four basic species groups: 
sandy faunal groups composed of species inhabiting both coarse and fine sand sub- 





Figure 2. Dendrogram resulting from normal classification analysis of 16 sampling local¬ 
ities at Barnstable Harbor. First characters indicate sampling station (see Fig. 1) ; last 
digit indicates sample number. 
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Figure 3. Dendrogram resulting from inverse clustering of species at the sampling stations. 
Species are indicated by code numbers (see Table II). 


strata; ubiquitous forms characterizing the dominant faunal components of Barn¬ 
stable Harbor; muddy and sandy-mud species; and single species groups which are 
both rare and restricted to a single sediment type. 

Animal-sediment relationships. While more detailed information will be pre¬ 
sented in a later paper (in preparation) on the feeding ecologies of the dominant 
benthic macrofauna, it is interesting to examine the general trophic-group sedi¬ 
ment-type relationships found at Barnstable Harbor. Table I lists information 
regarding the abundance and species number of the two major trophic groups: 
suspension- and deposit-feeders. The majority of the benthic fauna at Barnstable 
Harbor are comprised of deposit-feeding species. Mya arenaria and Gemma gemma 
are the two most abundant suspension-feeding organisms. These forms dominate 
at stations 12-15 which contain well-sorted, relatively fine sands with low silt- 
clay and organic carbon values. G. gemma comprises over 99% of the abundance 
of the suspension-feeding forms at these stations. G. gemma is widely distributed 
and abundant in fine sands. In coarser sand localities, where surface rippling 
becomes increasingly pronounced, the bivalve density is low. Sediment instability 
resulting from waves and currents readily move this small bivalve and are probably 
responsible in determining its distributional pattern. M. arenaria is much rarer 
and comprises less than 1% of the total number of suspension-feeders at the 
sampling stations. Bradley and Cooke (1959) and Sanders et al. (1962) have 
hypothesized that high density populations of G. gemma inhibit larval settling of M. 
arenaria. Deposit-feeders predominate at stations 1-9 where both organic matter 
and silt-clay fractions are the highest. The coarse sand, ripple-marked stations 
(10, 11, 16) show reduced numbers and species diversity of both trophic groups. 

Seasonal changes in community structure and species composition 

Station 1 was sampled at approximately monthly intervals for a period of 19 
months to obtain information on the degree of temporal fluctuation in both com¬ 
munity structure and species composition, 
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While the basic unit of structure in any natural community is the individual, 
examination of the total species complement of a particular assemblage or com¬ 
munity of organisms can provide insight into distinguishing temporal changes of 
community organization. Figure 4 summarizes three such measures of com¬ 
munity structure at station 1 : species diversity (H'), evenness (J'), and 
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Figure 4. Temporal changes in number of species/245 cm 2 (A) ; evenness (B) ; and 
species diversity (C) at station 1. 
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number of species, While all three measures show some temporal variation 
throughout the sampling period, they did not vary significantly according to the 
results of one-way analysis of variance (P > 0.05). It should he noted, however, 
that indices similar to the ones used are not sensitive to subtle changes in species 
composition. In particular, Shannon’s index H' is biased toward the more dom¬ 
inant species and the rarer members contribute little to the measure. Also, 
these measures do not provide any information regarding what species are actually 
characteristic of a given sample(s) during different periods of the year. For 
these reasons classification analysis was adopted to examine the degree of temporal 
partitioning of samples collected throughout the year and to elucidate changes in 
species composition at station 1. 

Normal analysis. Figure 5 presents the results of normal classification analysis 
of 46 samples collected from June 1974 to December 1975. Generally the major 
groupings derived from this procedure are subdivided by seasonal similarity and 
can he grouped into four clusters. Group A is largely a “winter’’ cluster comprising 
samples collected from October to January. Group P> contains collections from 
March-May and can he regarded as a “spring” cluster. Groups C and D can 
he classified as “summer” and “fall” groups, respectively, though the level of sepa¬ 
ration of these two clusters is higher than found in Groups A and B. 

It is evident that there is a certain degree of sample overlap between the clusters, 
indicating spatial variation in species abundance and composition of individual 
samples. Group B (spring cluster) contains one sample collected during the fall 
(October, 1975). Groups C and D have the greatest amount of sample overlap 
between clusters; each contain two samples that are not similar on a seasonal basis 
with other samples in the group. In all cases, upon examination of the original 
data matrix these samples were characterized by the absence of the rarer species 
which normally had occurred at other times of the year in the other samples. 


O = M-OoOWWW^^(D^- n !^^ > j00(DO03_® 
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Figure 5. Dendrogram resulting from normal analysis of sampling station 1 from June, 
1974 to December, 1975. First characters indicate the month sample was collected (prime 
numbers = 1975). Last number indicates sample number. 
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Inverse analysis. Species that occurred less than five times during the sampling 
program were excluded from the inverse analysis. Those species eliminated 
included six species that occurred only once, two forms occurring twice, and one 
species occurring three times. Twenty-one species were included in the analysis. 

Again, four basic seasonal species groupings are evident at station 1 (Fig. 6). 
As would be expected, individual clusters are composed of species whose seasonal 
peak population abundance levels coincide with one another. Table III provides 
information on the time of the highest numerical abundance coupled with data 
on the reproductive activities of each species included in the analysis. Data on 
reproductive behavior of each species were compiled from published sources and 
personal observations. Unfortunately, some of these data are based on limited 
observations. In general, however, there is a close concordance between the peak 
population abundance levels of species and their periods of reproductive activity. 
In most instances there is a one to two month delay between the observation of 
reproductively active adults and the appearance of the new year class into the 
benthic population. This time lag may either be associated with the planktonic 
larval phase exhibited by most of the species in this study or may be the result 
of sampling bias. Although the use of smaller meshed screens proved impractical 
in this study, they may provide a more exact time scale of larvae entering the 
benthic populations. 

Species group 1 is composed of organisms whose recruitment generally occurs 
in the spring and early summer months. This group contains those faunal ele¬ 
ments which are normally the dominant components of this station. Group 2 
contains species that are generally most abundant during the warmest fall months 
of the year. Four species of polychaetes and one species of amphipod, whose peak 
abundances occur in the coldest months of the year, comprise Group 3. Group 4 
contains six species whose recruitment of individuals into benthic populations 
occurs during the late summer months of July-September. 

Population charactcristcis of the dominant species. Figure 7 shows a sum¬ 
mary of the seasonal changes in the population abundance of five of the most 
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Figure 6. Dendrogram indicating seasonal species groups at station 1 using inverse 
classification analysis. Species are indicated by code numbers (see Table II). 
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numerically important species at station 1. S. bcncdicti, generally the most 
abundant and ubiquitous polychaete found in Barnstable Harbor, comprises 20- 
70% of the total fauna at station 1. Abundance patterns of this small spionid are 
highly seasonal with greatest numbers occurring in early spring. Rapid decreases 
occurred in the summer months in both 1974 and 1975. Much of this decline is 
attributable to the loss of larger individuals from the population. vSize-frequency 
distributions of this organism are uni-modal, characteristic of annual species. The 
new-year class enters the population in June and continues to grow throughout 
the winter and following spring (unpublished). Watling (1975) and Jones 
(1961) found similar seasonal abundance patterns of S. bcncdicti in Delaware 
and southern California populations, respectively. Watling (1975) listed size data 
indicating the annual nature of this species. 

E. hctcropoda and H. filiform is also appear to be relatively short-lived. Size- 
frequency distributions of E. hctcropoda indicate three separate modes, suggesting 
at least three year classes are found in the population. The 0-year class enters the 
population in July and continues to grow throughout the year. The population 
abundance of IE filiformis , on the other hand, was quite stable throughout the sam¬ 
pling program. Changes in the size-frequency distributions of this species indicate the 
0-year class enters the benthic population in July-August. The size distributions 
also indicate two distinct size classes, suggesting first and second year groups. 
Buchanan and Warwick (1974) found a similar life history pattern in IE fili - 
formis; however, 0-year groups appeared in January on the coast of Britain. 



1974 1975 

Figure 7. Estimated abundance (number per 245 cm 2 ) of the five most common macro- 
benthic species at station 1 from June 1974 to December, 1975. 
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Table III 


Seasonal species groups at Barnstable Harbor, time of peak population abundance, reproductive 
activity , and geographic range along the eastern coast of North America. 


Group 

Species 

Peak Abundance 

Reproductive activity* 

Range** 

1 

Streblospio benedicti 

J line 

June-Oct. 

c 


Heteromastus filiform is 

July-Aug. 

June-J uly 

c 


Scolelepis squamata 

April-June 

Mar.-J une 

c 


Lumbrinereis tenuis 

Alay-June 

June 

c 


Spiophanes bombyx 

May-June 

I une-Aug. 

c 


Eteone heteropoda 

J une—J uly 

June (?) 

c 


Gemma gemma 

May-June 

May-Aug. 

c 

2 

Tharyx sp. 

Sept.-Oct. 


V 


Nereis arenaceodonta 

Aug.-Oct. 

Aug. 

V 


Spio setosa 

Sept. 

Oct.-Feb. 

V 

3 

Scolecolepides viridis 

Feb.-May 

Oct.-Nov. 

N 


Scoloplos acutus 

Feb.-April 

Feb.-July 

N 


A ricidea jeffreysii 

April-May 

July 

C 


Paraonis fulgens 

Jan.-Feb. 


N 


Corophium insidiosum 

April-May 

April-Aug. 

N 

4 

My a arena ria 

Sept. Oct. 

Sept.-Oct. 

C 


Nassarius obsoletus 

Aug.-Sept. 

June-J uly 

C 


Retusa canalicidata 

\ug.-Sept. 

June-Aug. 

C 


Prionospio heterobra nchia 

July 

April—J uly (?) 

Y 


Capitella capitata 

July-Aug. 

April-Nov. 

C 


Pol in ices duplicatus 

Aug.-Sept. 


V 


* Information obtained from: personal observations; Simon and Brander, 1967; Pettibone, 
1963; George, 1966; Dean, 1965; Grassle and Grassle, 1974; Scheltema, 1962; Turner, 1951; 
Grant, 1965. 

** Biogeographic range of species: C = generally ranging along entire North American coast; 
Y = generally occurring from Massachusetts to Florida; N = normally ranging from Maine to 
Massachusetts. Information from: Pettibone, 1963; Day, 1973; Foster, 1971; Abbott, 1954; 
Bousheld, 1973. 

The large spionid, Spio sclosa, showed two separate population abundance 
peaks, one which occurred in November and a slightly larger one which occurred 
in May (Fig. 7). S. sclosa appears to exhibit two reproductive periods: first, 
a spring period in which eggs are laid in female tubes and develop into benthic 
larvae; and secondly, a fall period in which pelagic larvae are abundant (Simon, 
1967, 1968). The two peaks in abundance found at station 1 probably reflect 
the two recruitment periods. 

Figure 7 shows that Scolccolcpidcs viridis is most abundant during February— 
May. During the remainder of the year, a small and relatively stable population 
level is found at station 1. George (1966) reports that N. viridis spawns from 
October—November in Nova Scotia, Canada, waters, and reproductive activity 
appears to be inversely correlated with increases in water temperature. 

Discussion 

As with other benthic studies, sedimentary parameters appear to be influencing 
general trophic and species composition of the macrobenthos at Barnstable Harbor. 
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The work of Sanders (1958, 1960) has brought into focus the importance of sedi¬ 
mentary features that influence the distribution of suspension- and deposit-feeding 
benthos. From examination of the theoretical properties of sediment dynamics, 
Sanders postulated a value of 0.18 mm as the optimal grain size for supporting 
suspension-feeders. He also suggested that high silt-clay fractions in sediments 
(relating to greater amounts of food) helped to explain the presence of the 
deposit-feeding benthos. Though these predictions have been successfully tested 
in the majority of subtidal surveys (but see McNulty, Work and Moore, 1962), 
they show inconsistencies when applied to intertidal studies (Sanders, ct al., 
1%2; Bloom, ct al., 1972). Examination of Table I indicates that areas dominated 
by suspension-feeders (though they normally comprise 1-2 species) have an average 
median grain size of 0.28 mm, larger than Sanders’ predicted value. The impor¬ 
tance of high silt-clay fractions influencing the distribution and abundance of deposit- 
feeders, however, is very evident in this table. 

Sanders (1958) also suggested that the distribution of suspension-feeders is 
greatly effected by sediment instability (median grain size greater than 0.1S mm). 
Data from most coarse sand stations (10, 11, and 15) at Barnstable Harbor indicate 
that the abundance of both deposit- and suspension-feeders is sharply reduced. 
These stations are characterized by pronounced surficial ripple marks, evidence of 
sediment instability. This instability coupled with low amounts of organic carbon 
and silt-clay fractions explain the reduced abundances of both trophic groups in 
these areas. 

The survey of Sanders ct al. (1962) of five locations on or near Huekins Island, 
Barnstable Harbor, found that two sediment types (muds and stable sands) 
were dominated bv deposit-feeding trophic biomass at four of the stations and did 
not fit the predictions stated above. Re-examination of the data, however, indicates 
that the numerical abundance of the two trophic groups tends to be correlated 
with amounts of silt-clav. Stations A and C (Sanders ct al., 1962) dominated by 
G. gemma , have the lowest amount of silt-clay. Deposit-feeders dominate at those 
stations with higher percentages of silt-clav. Median grain size relationships, 
however, show no significant trends in predicting benthic trophic structure. The 
reasons for the discrepancies in the two sets of data are unknown. Comparisons 
of both sedimentary parameters and rank order abundance of the macrobenthos 
indicate major long-term (over 15 years) changes at all localities have taken place. 

While examination of temporal changes at one sedimentary type at Barnstable 
Harbor revealed distinct seasonal change in species composition, no patterns in 
species diversity, evenness, or number of species per sampling date occurred. 
Normal classification analysis indicated seasonal clustering of samples collected over 
a 19 month period. Clusters, though exhibiting overlapping elements between these 
groups, show repeatable yearly trends characteristic of habitats in a dynamic 
“equilibrium,” not successional state. The majority of benthic species at Barnstable 
Harbor exhibit distinctive seasonal fluctuations in numerical abundance, and most 
appear to be rather short-lived (1-3 years). Due to the dynamic temporal 
nature of recruitment and mortality patterns characteristic of short-lived species, 
it is not surprising that one finds seasonal species specific population trends at 
station 1. Inverse classification analysis clustered species into four groups, 
generally corresponding to seasonal population abundance peaks and periods of 
reproductive activity. 
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The majority of studies examining mechanisms influencing seasonal abundance 
patterns of species attribute such change to regular or systematic change in resource 
quality and abundance (Fretwell, 1972). Since the benthic fauna inhabiting the 
intertidal sand and mud flats are greatly affected by sedimentary characteristics, 
seasonal changes in sediment could promote population change. As mentioned, 
however, there were no significant changes in sediment size composition over a 
19 month period. Thorson (1957) and Minis (1967) have suggested that settling- 
patterns of many benthic organisms coincide with food availability. Seasonal 
fluctuations in planktonic organisms arc common in this area of New England 
(Lillick. 1937) and could affect benthic organisms feeding upon them. However, 
most of the species occurring at Barnstable Harbor are infaunal deposit-feeders 
which obtain the majority of their nutrition directly from sediments. Temporal 
change in sediment food supply (measured as total organic carbon and potential food 
particulate matter) remained relatively stable throughout the year at most of the 
sampling stations (Whitlatch, in preparation). The effects of well-buffered seasonal 
food supplies influencing population characteristics of deposit-feeders has been 
reviewed by Levinton (1972). 

Temperature is commonly considered to be the most important single factor 
influencing the distribution and reproductive activities of marine invertebrates 
(Thorson, 1946). Temperature fluctuations at Barnstable Harbor show marked 
temporal variation and are repeatable from year to year; two features that could 
affect observed repeatable seasonal clustering patterns of species groups. While 
temperature may provide a seasonal cue for community change, it does not totally 
explain why some groups of organisms are commonly found in the winter months 
while others are restricted to the summer months. 

A partial explanation of the occurrence of seasonal species groups can be ob¬ 
tained from a closer examination of the latitudinal distributional patterns of many 
of the species. Cape Cod, Massachusetts, has long been considered a zoogeographic 
boundary between the southern part of the cold-temperate Nova Scotian faunal 
province and the northern portion of the mild-temperate Virginian faunal province 
(see Hazel, 1970, for a historical review) due to differences in summer water 
temperature found on the north and south sides of the Cape (Schopf, 1967). 
Several studies of offshore benthic organisms have shown distinct faunal elements 
characteristic of these two areas (Parker, 1948; Peterson, 1964; Schopf, 1965; 
Hazel, 1970). Barnstable Harbor is somewhat unique in that it contains faunal 
elements common to both zoogeographic provinces, coupled with a group of species 
found to occur throughout the entire eastern coast of North America. 

Distributional ranges of the various species, summarized in Table III, indicate 
several distinct trends when related to seasonal temperature fluctuations. The 
Virginian components are generally found at their greatest abundance during the 
warmest months of the year (August-October), while the Nova Scotian elements 
are most common during the coldest temperature months of the year (January- 
April). The cosmopolitan forms, comprising the largest group of species, generally 
are most abundant during early spring and summer months. The Nova Scotian 
and Virginian faunal elements do not seem to be geographically limited by winter 
temperature since shallow water temperatures are about the same on both sides 
of the Cape (Schroeder, 1966; Schopf, 1967). Potentially southward migrating 
forms are limited by summer temperatures too high for adult survival. North- 
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ward migrating species are limited by summer temperatures that never are high 
enough for reproductive activity to he initiated. Summer water temperatures at 
Barnstable Harbor are higher than surrounding Cape Cod Bay water due to the 
shallow nature of the area, explaining the occurrence of Virginian faunal com¬ 
ponents. 

While the examination of faunal elements characteristic of the two zoogeographic 
provinces helped to shed light on mechanisms influencing some seasonal community 
change, other species may be responding to the effects of biological interactions 
(mainly competition and predation) from other species. Certain species of deposit¬ 
feeding polychaetes which overlap greatly in both food and space utilization are 
temporally segregated, thus reducing the possibility of competitive interactions 
in resource use (Whitlatch, in preparation). Also, common epifaunal predators 
at Barnstable Harbor (mainly Fitmlithts hctcroclitus and Nassarius obsoletus) 
which are abundant during the summer-fall months, can exert disproportionate sea¬ 
sonal effects on infaunal populations and can alter the species dominance struc¬ 
ture of soft-bottom communities (Schneider, 1976; Whitlatch, unpublished). 
Species reproducing during times of minimal predator activity could significantly 
reduce tbe effects of predation. 

The degree to which seasonal clustering patterns at the various habitats in 
Barnstable Harbor persist for longer periods of time is unknown at present. 
Comparisons of the present study with other studies of Barnstable Harbor 
(Sanders, ct al., 1962; Mills, 1967) revealed major long-term changes in both 
sedimentary features and species structure at several localities during the summer 
months. While temperature appears to be one of the controlling mechanisms for 
initiating reproductive activity in most species, the exact timing of the appearance 
of larvae in benthic populations certainly is influenced by year-to-year fluctuations 
in temperature cycles characteristic of this area of New England. This unpredict¬ 
able environmental variable will place constraints on the degree of long-term 
stability of the re-occurring seasonal species groups at Barnstable Harbor. Lie 
and Evans (1973) provide one of the few studies examining long-term (7 years) 
changes in species composition of various faunal assemblages in subtidal areas of 
Puget Sound, Washington. While they found relatively small species changes, 
Lie (1968) indicated that there was little fluctuation in environmental parameters 
( c.g. } temperature) in this area—a feature unlike Barnstable Harbor. 

Seasonal changes in community structure have been noted in several subtidal 
surveys ( e.g ., Tenore, 1972; Boesch, 1973) and several authors have recorded 
seasonal groups of particular organisms in subtidal areas (Levings, 1975; Muus, 
1967; Bodiou and Chardy, 1973 ; Grassle and Smith, 1976). Environments char¬ 
acteristic of seasonal patterns of species composition are relatively unstable and 
“physically-controlled” (Sanders, 1968) ; namely, estuarine or intertidal and shal- 
low-subtidal habitats where both the probability of environmental disturbance and 
magnitude of seasonal fluctuations in environmental parameters are more 
likely to affect biological changes. The distinct seasonal clusterings of species at 
Barnstable Harbor, while an important characteristic of the area, may be even 
more pronounced than other areas owing to its somewhat unique geographic posi¬ 
tion and hydrographic conditions which have the potential of magnifying the 
effects of temporal community change. 
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In a review of tlie community concept and its validity when applied to marine 
benthic ecology, Stephenson (1973) pointed out the lack of the addition of a 
temporal component to most community definitions. While most benthic ecologists 
view communities as loosely organized suites of rc-occurring species, they fail to 
recognize the dynamic temporal nature of the organization of the component 
members. This failure can be largely attributed to the short-term (“one-shot”) 
nature of sampling in the vast majority of sampling programs. To better under¬ 
stand the nature of communities, physical and biological parameters which in¬ 
fluence them, and their persistence or stability, it is important to gain knowledge 
of temporal changes in community structure. Understanding these patterns will 
provide further insight into the nature of the structure and evolution of marine 
soft-bottom biological systems. 

This paper forms a portion of a thesis submitted to the University of Chicago 
in partial fulfillment of the requirements for a Ph.D. degree. I would like to 
dedicate this work to Ralph G. Johnson, friend and advisor, whose significant con¬ 
tributions in the field of marine benthic ecology have been a constant source of 
personal stimulation. Dr. IT. L. Sanders and J. F. Grassle provided important 
insight into various aspects of the study. Dr. Grassle provided use of the classifica¬ 
tion program. I wish to thank the assistance of Ty Ranta, Barnstable Conserva¬ 
tion Department for unlimited access to the marsh and use of a boat for collecting 
samples. L. F. Boyer, G. Johnson, and A. Williams aided with field work. Grants 
from Sigma Xi and the Hinds Fund for Evolutionary Research, University of 
Chicago, helped in starting the project. This research was largely supported by 
nsf grant I)ES 72-0160S-A01 and an equipment grant from the Louis B. Block 
Fund made to R. G. Johnson, University of Chicago. 

Summary 

1. A quantitative sampling survey of the benthic macrofauna inhabiting the 
intertidal sand and mud flats of Barnstable Harbor, Massachusetts, was conducted 
to describe general community structure and examine temporal changes in species 
composition. 

2. Classification analysis delimited coarse and fine sand, mud, mtiddy-sand, 
and gravel-mud benthic species associations. The 32 species used in the inverse 
classification analysis were partitioned into 10 species groups, reflecting spatial 
distributional patterns. Many of the species were both dominant and ubiquitous, 
masking discrete species groupings. 

3. The majority of macrobenthos at Barnstable Harbor were deposit-feeders 
which comprised more than 90G of all organisms sampled. The deposit-feeders 
normally dominate mud and muddv-sand sediments. Suspension-feeders were 
most abundant in fine sands. The relationship of sedimentary parameters affecting 
the distribution of both trophic groups proposed by Sanders is generally supported. 

4. While no significant changes were evident in species diversity, evenness, or 
species number throughout a 19 month sampling period, classification analysis de¬ 
limited seasonal clustering of both samples and species groupings. These pat¬ 
terns were repeatable over a two year period suggesting that dynamic “equilibrium,” 
not successional change was occurring. 
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5. Seasonal clustering* patterns of some species were related to the appearance 

of different zoogeographic faunal province elements. Typical warm-temperature 

Virginian components and cold-temperate Nova Scotian species were commonly 

found in late summer and winter, respectively. 
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